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ABSTRACT 



> 

^sq ■ We describe a possible new technique for precise wavelength calibration of high-resolution 

^vq I astronomical spectrographs using femtosecond-pulsed mode-locked lasers controlled by sta- 

\^ . ble oscillators such as atomic clocks. Such 'frequency combs' provide a series of narrow 

CO ' modes which are uniformly spaced according to the laser's pulse repetition rate and whose 

\ absolute frequencies are known a priori with relative precision better than 10"^^. Simula- 

te ■ tions of frequency comb spectra show that the photon-limited wavelength calibration preci- 

' sion achievable with existing echelle spectrographs should be ~ 1 cm s"^ when integrated over 

a 4000 A range. Moreover, comb spectra may be used to accurately characterise distortions 
Qh" of the wavelength scale introduced by the spectrograph and detector system. The simulations 

JL I show that frequency combs with pulse repetition rates of 5-30 GHz are required, given the 

^ ■ typical resolving power of existing and possible future echelle spectrographs. Achieving such 

^ ■ high repetition rates, together with the desire to produce all comb modes with uniform inten- 

^ ! sity over the entire optical range, represent the only significant challenges in the design of 

a practical system. Frequency comb systems may remove wavelength calibration uncertain- 
ties from all practical spectroscopic experiments, even those combining data from diff'erent 
telescopes over many decades. 



>< 



^ ■ Key words: instrumentation: spectrographs - instrumentation: detectors - methods: labora- 

tory - techniques: spectroscopic 



1 INTRODUCTION 

Echelle spectrographs are the basic tool for astronomical spec- 
troscopy with high velocity precision. Since many echelle diffrac- 
tion orders can be cross-dispersed and thereby recorded simulta- 
neously on rectangular-format media [e.g. modern charge-coupled 
devices (CCDs)], more-or-less continuous wavelength coverage 
over much of the optical range can be achieved whilst maintaining 
high resolving powefl Most modern echelle spectrographs operate 
at = 30000-150000; higher resolving powers usually entail re- 
duced wavelength coverage. However, to fully exploit the velocity 
precision available in high-resolution spectra, precise wavelength 
calibration over much of the available wavelength range is usually 
required. For example, se arches for extra- solar p l anets via the ra- 
dial velocity method (e.g. iMavor & Oueloz|[l995l : iMarcv & Butlej 
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^ Resolving power R = /1/FWHM, where FWHM is the instrumental pro- 
file's full-width at half-maximum. 



[1993) have achieved night-to-night relative precisions better than 
lms~^ from stellar absorption lines integrated over wavelength 
ranges of ~1000-3000-A (e.g. iRupprecht et al.ll2004l : iLovis et all 
2006). Also, echelle spectra of quasar absorption systems have con- 
strained cosmological variations in fundamental constants, such 
as the fine- struct ure constant (e.g. iBahcalL Sargent & Schmidtl 
19671: IWebb et al. 1999) and pro t on-to-electron mass r atio (e.g. 



Varshalovich & Levshakovl 1 19931 : llvanchik et all l2005h . Recent 



measurements of of-variation using many absorption systems 
reach relative precisions of ~10~^ when comparing severa l tran- 
sitions separated by up to -500 A (e.g. Murphy et al] l200ll : 
iMurphv. Webb & Flambauml l2003h . The corresponding velocity 
precision required is ~20 m s~^ 

The need for even higher precision is already clear. Detect- 
ing Earth-mass extra-solar planets around solar-mass stars will re- 
quire <10cms~^ precision stable over many-y ear timescales . Also , 
current evidence for varying constants (e.g. iMurphv et al.ll2004l : 
iReinhold et al.l [200^ remains controversial and such potentially 
fundamental results should be refuted or confirmed with the high- 
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Table 1. The requirements for an ideal wavelength calibration source and 
the relative advantages and disadvantages of three current sources plus the 
frequency comb possibility we describe in this paper. We consider a spectro- 
graph with resolving power R ~ 150000. 'Poss.' means possible depending 
on the design and construction of the calibration source. 



Calibration lines 


ThAr 


I2 cell 


Etalon 


Comb 


From fundamental physics 


Yes 


Yes 


No 


Yes 


Individually unresolved 


Mostly 


Yes 


Poss. 


Yes 


Resolved from each other 


No 


No 


Poss. 


Poss. 


Uniformly spaced 


No 


No 


Yes 


Yes 


Cover optical range 


Yes 


No 


No 


Poss. 


Uniform strength 


No 


No 


Poss. 


Poss. 


Long-term stability 


No 


Poss. 


No 


Yes 


Maintain object S IN 


Yes 


No 


Yes 


Yes 


Exchangeable 


Yes 


Yes 


Poss. 


Yes 


Easy to use 


Yes 


Yes 


Poss. 


Poss. 


Reasonably low cost 


Yes 


Yes 


Yes 


Poss. 



est possible confidence. There is even the possibility of measur- 
ing the acceleration of the Universal expansion by monitoring the 
change in r edshift of Lyman-cy forest a bsorption lines in quasar 
spectra (e.g. ' Sandagelll962l : lLoeblll998l ). This would require two 
epochs of spectra, taken decades apart, with velocity precision ap- 
proaching 1 cms~^ integrated over ~ 1000-3000 A ranges. This is 
one possible aim of the Cosmic Dynamics Experiment (CODEX) 
spectr ograph proposed for the European Extre mely Large Tele- 
scope /Pasauini et al.l l2006l : iGrazian et alJl2007l : Liske et al. 2007, 
in preparation). 

Astronomical echelle spectra are usually wavelength cali- 
brated using comparison emission- or absorption-line spectra. The 
ideal comparison spectrum would comprise a series of lines which: 
(i) have known wavelengths which are determined by fundamen- 
tal physics; (ii) are individually unresolved; (iii) are resolved from 
each other; (iv) have uniform spacing; (v) cover the entire opti- 
cal wavelength range; (vi) have uniform strength (intensity or op- 
tical depth); (vii) are stable over long time-scales and (viii) do 
not reduce the signal-to-noise ratio {S IN) of the object spectrum. 
The calibration source should also be: (ix) exchangeable, in the 
sense that two independent sources produce the same spectrum; 
(x) easy to use, i.e. essentially 'turn-key' operation for non-expert 
users and (xi) should have a reasonably low cost. Table [T] sum- 
marises whether four different calibration sources meet these re- 
quirements. The most commonly used calibration source is the 
thorium-argon (ThAr) hollow-cathode emission-line lamp. While it 
is relatively low cost, has lines whose wavelengths are determined 
by atomic physics and which cover the entire optical wavelength 
range, it has several disadvantages. For example, ThAr lines differ 
widely in intensity and spacing. Indeed, most are blended together 
dX R < 150000. Furthermore, as the current and pressure in the 
lamp varies (e.g. as the lamp ages) the relative intensity of differ- 
ent ThAr lines varies widely. Iodine cells are also commonly used 
in radial-velocity extra-solar planet searches (e.g. IMarcv & Butlej 
11992). These have similar advantages and disadvantages as ThAr 
lamps but they are typically only useful over the wavelength range 
5000-6500 A and absorb about half the object light, thereby signif- 
icantly reducing the S IN. 

A significant step towards the ideal calibration source might 
be achieved with the relatively new technology of optical 'fre- 
quency combs' generated from mode-locked femtosecond-pulsed 
lasers. The 2005 Nobel Prize in physics was awarded for the pio- 



neering work on such combs (e. g. 'Reiche rt et aLlI 19991 : Ijones et aP 
I2OOQI : lUdem. Holzwarth & Han sch 2002). As the name suggests, 
a frequency comb provides a spectrum of continuous-wave laser 
lines called modes. They are produced by the repetitive pulse train 
of the mode-locked laser. The mode spacing, which is constant in 
frequency space, is given by the pulse repetition frequency. Since it 
resides in the radio frequency domain it can readily be synchronised 
with an extremely precise frequency reference such as an atomic 
clock. The comb's perfectly regular frequency grid should also fa- 
cilitate the effective compensation of non-linearities in the spectro- 
graph's calibration curve while directly linking the calibration to 
the fundamental physics of the atomic clock. This direct link to the 
definition of time and the long-term stability mean that even if dif- 
ferent combs are used to calibrate spectra taken at (perhaps widely) 
different epochs, those spectra can be directly compared. Research 
utilising archival and 'virtual observatory' data would directly ben- 
efit from such a reproducible calibration source. These properties 
would greatly facilitate extra- solar planet detection and varying- 
constant analyses and would be more-or-less required for detect- 
ing the Universal acceleration via Lyman-of forest spectra taken 
decades apart. 

This paper describes further the possibility of frequency comb 
calibration of echelle spectrographs. Section [21 outlines the basics 
of frequency comb generation, the methods used to ensure that they 
are stable and the general manner in which they might be used as 
practical echelle calibration sources. Section [3] derives the optimal 
laser repetition frequency and the overall precision of the wave- 
length calibration that might be achieved with typical echelle spec- 
trographs. Some likely sources of systematic errors which may oth- 
erwise limit the precision are also discussed. Section|4l outlines the 
main challenges in designing a working frequency comb calibra- 
tion system and proposes some possible solutions. 



2 BASIC OPERATION 
2.1 Generating the comb 

The exact properties of a frequency comb may be derived even 
without knowledge of the operational principles of a mode-locked 
laser. It suffices to know that these lasers can store a single pulse 
and maintain it on a repetitive path. After each round-trip a copy of 
this pulse is emitted though the output mirror/coupler of the laser, 
resulting in an indefinite train of laser pulses. Like in any other 
laser, the energy lost is replenished by stimulated emission in the 
lasing medium, thereby allowing the energy in a short pulse to be 
stored for any length of time by compensation of loss with gain. 

To understand the mode structure of a femtosecond frequency 
comb and the techniques employed to stabilize it, consider the ide- 
alized case of a pulse circulating as a carrier wave with frequency 
Vc in a laser cavity with length L . The cavity strongly modulates 
the carrier wave amplitude which can therefore be described by a 
periodic envelope function A(t) = A(t - T), where T is the round- 
trip travel time of the pulse calculated from the cavity mean group 
velocity: T = ILIVg^- The repetition rate of the laser is simply 
Vr = 1/r. This is illustrated in Fig. [T] Since the envelope func- 
tion is periodic it can be written as a Fourier series and the electric 
field at any given place (e.g. the laser's output mirror/coupler) can 
therefore be written as 

E(t) = A(t)e-'^'''^' = Am^"''^^^"'^^"^^' , (1) 

m 

where A,„ is the mth Fourier component of A{t). Equation ([T]) shows 
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Figure 1. A pulse train (top) produces the optical frequency comb in Fourier 
space (bottom) with some properties inverted. Shorter pulse envelopes pro- 
duce broader frequency combs. The round-trip time of the pulses inside the 
generating laser cavity (not shown) determines the repetition frequency, Vr, 
exactly. Dispersive elements in the laser cause a difference in the group and 
phase velocity, continuously shifting the carrier wave with respect to the 
envelope in the upper trace by per pulse. In the frequency domain this 
causes the comb to shift by an amount Vce = iS^ipjlnT . 



that, under the assumption of a periodic pulse envelope, the result- 
ing spectrum is a comb of laser frequencies with spacing Vr such 
that the mth frequency component is = mvr -h Vc. In general, Vc 
is not an integer multiple of Vr and so the modes are shifted from 
being exact harmonics of Vr. Re-numbering the modes gives 



(2) 



for n SL large integer, where the carrier envelope offset frequency, 
Vce, is restricted to < Vce < Vr- This equation maps two radio fre- 
quencies, Vr and Vce - which can be stabilized by a precise reference 
such as an atomic clock - onto the optical frequencies Vn. 

Given the simple treatment above, one might expect deviations 
from equation ^ in reality. However, mode locking the laser en- 
sures that equation (O is realised to very high precision in practice. 
The frequencies y„ in equation ^ represent the longitudinal cav- 
ity modes of the mode-locked laser. The process of mode locking 
stabilizes the mode frequency difference, Vr, and the cavity mode 
phases such that their superposition yields a short pulse. The exis- 
tence of such a pulse depends on the stability of the relative phases 
of the modes: if these change in time - for example, due to a small 
deviation of the mutual mode separation from - the pulse would 
disintegrate immediately. The very observation that the pulse stays 
intact for a time At implies that variations in Vr must be < For 
a pulse life-time of 1 s, this is only a few parts in 10^"^ at the optical 
frequencies Of course, the stored pulse lasts much longer than 
1 s. In fact, deviations from equation ^ - that is, deviations from 
a constant m ode spacing - have so f ar evaded detection even at the 
-10-16 level (iHolzwarth et"aDl200Ql) . 

As the pulse circulates in the cavity, dispersion causes a differ- 
ence between the group and phase velocities. The pulse envelope, 
A(0, propagates at the group velocity, Vgr, while the carrier wave 
travels at its phase velocity, causing them to shift with respect to 
each other after each pulse round-trip. This is illustrated in Fig.[T] 
Comparing two consecutive pulses with A(t) = A(t - T), equation 
([T]) implies that the electric field is given by 



E(t) = E(t - T)e- 



E(t - T)e- 



(3) 
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where, in the last step, a large integer multiple of In was subtracted 
from the phase, equation (O was used with Vc as one of the modes 



Figure 2. A typical self-referencing scheme. Top: The carrier envelope off- 
set frequency, Vce, that displaces the modes from being exact harmonics of 
the repetition rate, Vr, is measured by frequency doubling some modes at 
the 'red' side of the comb and beating them against modes at the 'blue' 
side (depending on the type of laser, these do not need to be truly blue or 
red). Bottom: Such a self referencing scheme employs a dichroic mirror to 
separate the 'red' and 'blue' wings of the frequency comb. The former is 
frequency-doubled in a non-linear crystal and reunited with the 'blue' part 
to create a wealth of beat-frequencies, all at Vce- 



and we recall that Vr = . That is, the phase shift for each pulse 
round- trip is = InVc^T . Therefore, the frequency comb is offset 
from being exact integer harmonics of Vr b y the so-called carr ier- 
envelope offset frequency, Vce = AipjlnT dUdem et all 12002b . It 
should be emphasised that the model for the origin of Vce is not 
important for precision measurements using the comb; all that is 
required is that Vce can be fixed to a certain value experimentally 
(e.g. via a feedback loop) and that the mode spacing, Vr, is constant. 

As Fig. [1] illustrates, the spectral envelope of the frequency 
comb has a width that is roughly given by the inverse of the pulse 
duration and a mode spacing given exactly by the inverse round- 
trip time, Vr = . Modern mode-locked lasers readily achieve 
~ 100 THz- wide spectra (=^830 A wide at 5000 A) from -lOfs 
pulses. Two types of lasers are common: the titanium-sapphire 
Kerr-lens mode-locked laser which operates around 8000 A and the 
erbium- and ytterbium-doped fibre lasers which operate at around 
1.55 and 1.04 |j,m respectively. Typical values for the mode spac- 
ing/repetition frequency are Vr ~ 100 MHz, but in some special 
cases up to Vr ~ 40-160 GHz has been achieved. See further dis- 
cussion in Section |4] 



2.2 Self-referencing the comb 

Measuring and/or stabilizing both radio frequencies, Vce and Vr, 
fixes the entire optical frequency comb by virtue of equation (|2]). 
While the pulse repetition rate, Vr, is readily measured with a photo- 
detector somewhere in the beam, the determination of Vce requires 
significant effort. In fact it took more than 20 years to develop a 
simple Vce detection scheme, dubbed 'self-referencing', which is il- 
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lustrated in Fig. [2 A low frequency mode with mode number ui is 
frequency-doubled to frequency 2{niVy- + Vce). This is superimposed 
with a high frequency mode, n2, onto a photo-detector. The result- 
ing beat-frequency is therefore (2«i - n2)Vx + Vce which reduces to 
the desired signal, Vce, provided that 2n\ - n2. This latter condition 
requires the comb to span an entire optical octave. Lasers capable 
of p roducing such broad spectra are only now beco ming available 
(e.g.lMatos et al. 2 004lFortier. Bartels & Diddams 2006): they are 
currently experimental and have several disadvantages. The stan- 
dard method of producing this spectral width has therefore been 
via a process called self-phase modulation induced by non-linear 
interactions in optical fibre^. For practical pulse repetition rates in 
the radio frequency domain, especially those approaching the GHz 
regime, the required peak intensity for such non-linear processes is 
difficult to achieve because the available laser power is shared by 
more pulses. However, with the development of Kerr-lens mode- 
locked lasers and photonic crystal fibres (e.g. Knight et al. 1996), 
octave- spanning combs can now be simply generated. The fibre 
lasers mentioned above use specially doped fibres to achieve simi- 
lar spectral broadening. 

Self-referencing applies the full accuracy of the radio fre- 
quency reference to the entire optical comb, yielding a perfectly 
regular calibration grid for the spectrograph. Current commercial 
caesium atomic clocks would provide a calibration as precise as 
several parts in 10^^, i.e. a velocity precision of ~0.01 cms~^ This 
may even be improved by using a timing signal broadcasted by the 
Global Positioning System to further stabilize the caesiuni cloc k 
or by setting up a local caesium fountain clock (e.g. lBauchll2003l) . 
However, we demonstrate in Section [3] that typical echelle spectro- 
graphs will limit the available precision to ~1 cms~^ and so such 
improvements should not be necessary. 

In contrast to this standard radio-frequency self-referencing, 
an alternative optical referencing scheme might prove more practi- 
cal for astronomical applications. In such a scheme the nth optical 
mode, y„, could be stabilized, for example, to an iodine- stabilized 
laser. The repetition rate would be stabilized to a precise radio fre- 
quency reference just as in the standard scheme above. As iodine- 
stabilized lasers can be as accurate as a few parts in 10^^ (e.g. 
lYe. Ma & Hal]l200lh this method could also provide the required 
accuracy. 

2.3 The comb as a calibration source 

Figure[2l shows the self-referenced (i.e. frequency calibrated) comb 
spectrum being fed into the echelle spectrograph via an optical fi- 
bre. We envisage that such a configuration would off'er the greatest 
stability and potential for tracking real-time distortions of the spec- 
trograph wavelength scale. In practice, the spectrograph would be 
fed by two fibres, one carrying the science object light and the other 
carrying the comb light for calibration. The two spectra would be 
recorded simultaneously next to each other on the CCD. The in- 
tensity of the frequency comb spectrum might be varied according 
to changes in the intensity of object light being transmitted by the 

^ Self-phase modulation originates from an intensity-dependent refractive 
index of the material. In particular, the polarizability of the material behaves 
non-linearly with respect to the incident electromagnetic field. This is also 
called the 'Kerr effect'. It provides an additional phase delay proportional 
to the field's intensity distribution (in both space and time). Applications in 
the context of the calibration comb include intensity-dependent lensing for 
mode-locking ('Kerr-lens mode-locking') in femtosecond lasers and spec- 
tral broadening or ultra-fast pulse shaping. 



telescope and spectrograph. Such variations could easily be tracked 
via feedback from detection of stray object light as the main beam 
passes into the spectrograph. Thus, changes in the weather con- 
ditions or small drifts in the spectrograph itself during the object 
exposure, which might ordinarily lead to distortions of the wave- 
length scale, would aff'ect the recorded comb spectrum in the same 
way. 

Of course, for very high precision work, one would have to 
consider other instrumental eff'ects, such as diff'erences between 
(and inhomogeneities in) the object and comb beam profiles as they 
emerge from the fibres. One would also consider placing the entire 
spectrograph in vacuum to ensure temperature and pressure stabil- 
ity. M any of the above techniques have already been demonstrated 
(e.g. Baranne et al.ll 19961 : iMavor et al.ll2003) . Of course, the more 
traditional approach of taking calibration exposures before and/or 
after the object exposure would still be open with the comb as a 
calibration source. The comb light could also be expanded and dif- 
fused after emerging from the fibre to illuminate a traditional spec- 
trograph slit. 

Note that the specifics of a final frequency comb design and 
spectrograph feed system need careful consideration quite beyond 
the scope of this paper. Indeed, these design considerations first 
require knowledge of the optimal laser repetition rate for typical 
echelle resolutions. The various tolerances involved in the design 
should also be informed by the photon-limited calibration preci- 
sion achievable with an echelle spectrograph. These quantities are 
calculated in Section[3]and their implications for design challenges 
are considered in Section |4l 



3 OPTIMAL PHOTON-LIMITED WAVELENGTH 
CALIBRATION PRECISION 

The wavelength calibration precision available from a comb spec- 
trum will clearly improve with increasing line-density. However, 
since the spectrograph resolving power is finite, too high a line- 
density will decrease the contrast between neighbouring lines and 
degrade the available precision; an optimal line-spacing must exist. 
Here we calculate this optimum laser repetition rate and the asso- 
ciated wavelength calibration precision. We first present a simple 
estimate of the photon-limited precision using some general argu- 
ments in Section 13.11 Sections 13.21 & 13.31 describe more rigorous 
calculations which also provide the optimal laser repetition rate. 
These calculations concentrate only on the photon-limited preci- 
sion and ignore systematic errors; the latter are discussed briefly 
in Section [331 We assume throughout this section that only one 
CCD exposure of the comb spectrum is used for wavelength cal- 
ibrating the spectrograph. The dynamic range of the CCD pixels 
therefore limits the comb S /N. Clearly, future improvements in de- 
tector technology, such as multiple or continuous reading out of the 
comb spectrum during a single science exposure, could improve 
the efl'ective comb S /N, thereby allowing even more precise cali- 
bration. 



3.1 General considerations 

Some general arguments serve as an intuitive introduction to the 
calculations in the following sections. Let us first consider a single 
unresolved comb line. When extracted from the CCD, this line will 
reflect the instrumental profile of the spectrograph. Assuming this 
to be a Gaussian whose FWHM is sampled by ;i > 1 pixels, a 
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well-known approximation to the error on the line's position is (e.g. 
iBrauldl 19871) 



dv^A 



FWHM 



(4) 



The S/N here refers to the peak S/N per pixel across the line 
(i.e. usually at the line centre). The pre-factor, A, depends on both 
the functional form of the line-profile and the relationship between 
S IN and pixel intensity. For example, if the noise level is constant 
across the pr ofile then A - 0.693 for a Gaussian and 0.798 for a 
Lorentzian ( Braultl[l987t) . In our case the noise is Poissonian with 
an additional small contribution from detector noise. From a sim- 
ple numerical experiment in which Gaussian profiles with this noise 
profile were fitted with standard least-squares techniques, we find 
that A ^ 0.41 when the photon noise dominates the detector noise. 

To determine the total velocity precision available from, say, 
an entire echelle order containing such comb lines we must know 
the line density and the spectral width of the order. As a first guess, 
let us assume that a comb line occurs every 2.5 resolution ele- 
ments so that the line-density is high but neighbouring lines do 
not strongly overlap. For consistency with the more rigorous cal- 
culations below, let us assume that echelle orders are 2048 pix- 
els long with 3-pixel sampling of the FWHM. That is, there are 
N ^ 272 comb lines per echelle order. Therefore, assuming a peak 
S IN - 500 in each line, the total calibration precision per echelle 
order at = 150000 (i.e . FWHM = 2kms-^) is a, ^ 6vl Vn = 
0.41 X 2kms-V(500 x V3 x 272) = 5.7 cm s-^ 

The following sections demonstrate that our assumption that 
comb lines should be spaced by ~2.5 resolution elements does in 
fact provide an optimal calibration with precision close to 6 cm s~^ 
per echelle order as predicted above. 

3.2 Velocity information content of comb spectra 

If a smal0 velocity change is to be measured from pixel / in a spec- 
trum F(i) w ith l-cr error array cr^ (/), the lim iting velocity precision 
is given by (iBouchv. Pepe & Ouelozll200"lh 



o-pii) 



(5) 



c A(i) [dF(i)ldA(i)] 

This equation simply states that a more precise velocity measure- 
ment is available from those pixels where the flux has a large gradi- 
ent and/or the uncertainty in the flux is small. We have used wave- 
length as the dispersion coordinate here but one may replace A with 
any other dispersion measure desired, including pixels. This quan- 
tity can be used as an optimal weight. 



W(i) = [o-,(i)lc]- 



(6) 



to derive the total velocity precision available from all pixels in a 
spectrum, 

' (7) 



If the spectram is divided into smaller spectral slices k, such as 
echelle orders, then it follows that the total velocity precision avail- 
able is a weighted average of that available from each slice. 



1 



(8) 



^ In all applications discussed in this paper, the expected velocity change 
is much smaller than the width of any spectral feature. 



3.3 Simulations of comb echelle spectra 

Equations & ([S]) can be used to calculate the velocity precision 
achievable with a frequency comb recorded with an echelle spectro- 
graph. To this end we construct a simple model of a frequency comb 
spectrum by writing the flux in each pixel, measured in photo- 
electrons, as the product of a slowly varying 'envelope function' 
Fe(/) and a quickly varying 'comb function' /(/) which varies be- 
tween zero and unity: F(i) = /(0^e(0- A reasonable model of the 
expected l-cr error array is then 



(9) 



where o-^{i) is the noise contribution from the detector in each pixel. 
The signal-to-noise ratio for each pixel can then be written as 



) = m 



m+Plii)' 



(10) 



where = o-d(i)l ^^FJJ) is the ratio of the detector noise and 
photon noise expected from the envelope at pixel /. Assuming that 
o'dii) = CTd is constant for all pixels and recalling that Fq(i) varies 
slowly compared to /(/) then the maximum S IN of the comb is 



(-) = 

VA^/max 



1 +/5d,max 



(11) 



Therefore, in our simulations one specifies the signal-to-noise at all 
pixels by knowing the shape of the comb and envelope functions 
and by specifying the maximum signal-to-noise, (*S/A/^)max, across 
the comb and a constant detector noise, cTd or /?ci,max- 

For our simulations we assume that the envelope function is 
constant over the entire wavelength range of the comb. That is, all 
comb lines have equal intensity. The comb function can be mod- 
elled as a series of delta functions equally spaced in frequency ac- 
cording to the assumed laser repetition rate, v^. The delta-functions 
are convolved with the spectrograph's instrumental profile which 
we assume to be a Gaussian of width specified by the resolving 
power R. For an echelle spectrograph R is roughly constant over 
the entire wavelength range and we make this additional assump- 
tion here. Simulated spectra are placed on a log-linear wavelength 
scale sampled at 3 pixels per resolution element (FWHM) in order 
to emulate one-dimensional reduced echelle spectra. A portion of 
simulated spectrum is shown in Fig. [3] 

We calculate o-^ from the simulated spectra using equation Q. 
In calculating the derivative in equation ^ we have simply used the 
finite-diff'erence derivative between successive pixels. 



dF(i) ^ F(i) - F(i - 1) 
dA(i) ^(0 - A(i - 1) 



(12) 



With this application, the first pixel of any spectral subdivision has 
undefined weight, W{i). The sum in equation Q therefore begins 
at the second pixel in each spectral subdivision. 

An example of the results from the simulated comb spectra 
is shown in Fig. U Here we assumed R = 150000 and (*S/A/^)max 
of 500 per spectral pixel. Gemini/bHROS, ESO-3.6m/HARPS and 
Subaru/HDS already operate at (or close to) this resolving power. 
To obtain (5'/A/^)max = 500 per pixel with modern CCDs, which 
typically have a linear response when fewer than ~40000 photo- 
electrons are detected per pixel, the comb spectrum would have to 
be spread over 6-7 pixels in the spatial direction; this is readily 
achievable with most modern echelle spectrographs. We also as- 
sumed a fairly typical wavelength range, 3800-8200 A. 
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Figure 3. Comparison of small sections of simulated object spectra with 
different calibration spectra. The upper panel shows a model Solar spec- 
trum (red/grey line) and a low-ionization metal transition (arbitrary red- 
shift) with typical velocity structure in a synthetic quasar absorption spec- 
trum (black line). The next panel compares model spectra of an iodine ab- 
sorption cell (red/grey line) and ThAr emission-line lamp (black line) at 
R = 150000 on different arbitrary flux scales. Note the complexity of the 
I2 absorption spectrum and the relatively sparse distribution of ThAr lines 
with their large range of intensities. The third panel shows a simulated fre- 
quency comb echelle spectrum over the same wavelength region with repe- 
tition rate Vr = 10 GHz, resolving power R = 150000 and (S /N)max = 500. 
The lower panel zooms in on a 1-A wide section of the upper panel. For 
the comb spectra, the black histogram shows the expected flux while the 
red/grey solid line shows the error array exaggerated by a factor of 25. 



Fig. m demonstrates two important results: (i) The velocity 
precision available from each echelle order is ~6-7cms~^ which 
is consistent with the 'back-of-the-envelope' calculation in Sec- 
tion [TT] This varies by <20per cent across all echelle orders be- 
cause the comb lines are uniformly separated in frequency space; 
(ii) The total velocity precision integrated over all echelle orders is 
0.45 cms~^ but, for a given R, this depends strongly on the comb 
line-spacing. Assuming that each spectral resolution element is 
sampled at or above the Nyquist rate and the detector noise is small 
compared to the photon noise, the optimum line- spacing required 
is determined by the resolving power: 



v7 



1.5 X 10^ 
R 



GHz ~ SAVce 



(13) 



where Avcent is the spectrograph's resolution in frequency space at 
approximately the centre of its wavelength range. Also, at constant 
R and (S /N)rriax (and again assuming /?ci,max ^1), the increase in 
overall velocity uncertainty with increasing is simply described 
as cTy oc provided that comb lines are well resolved from 
each other, i.e. Vr > 1.5yr^\ The increase in cry with decreasing 
line- spacing below v"^^ is extremely fast but has no simple ana- 
lytic description. Suffice it to say that our simulations show that the 
increase is faster than exponential. Nevertheless, a total precision 
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Figure 4. Results from comb simulations with R = 150000, (S/N)raax = 
500 and ySd,max = 0.01. The spectrograph's resolution (FWHM) is sampled 
with 3 pixels and each echelle order is assumed to be 2048 pixels long. 
The main panel shows the expected velocity precision available from each 
echelle order, o-y(k), for three difl'erent line-spacings while the inset shows 
the total velocity precision, cry, integrated over the full wavelength range of 
the simulated spectrum, 3800-8200 A. For a comb spectrum of this resolu- 
tion the optimum line-spacing (or repetition rate) is Vr = 10 GHz, resulting 
in a total velocity uncertainty of o-y = 0.45 cm s~^ o-y rises sharply for 
repetition rates below the spectrograph's resolution (3.3 GHz at 6000 A). 



within a factor of 2 of the optimum can be achieved over the range 
Vr ^ (0.5-5) v^^\ Equivalently, a comb of fixed Vr can be used to 
almost-optimally calibrate spectra over a reasonable range in re- 
solving power. 

These results are easily scaled according to the parameters 
of diff'erent spectrographs. It is not surprising to find that CTy oc 
\/(S /A/^)max and o-y oc \/Rif, again, the sampling of each resolution 
element is adequate and /?d,max ^ 1- Thus, for a given resolving 
power and maximum S /N, the optimum photon-limited precision 
available over the wavelength range 3800-8200 A is 



opt 



0.45 



0.45 



500 



500 



1.5 X 10^ 



R 

1.5 X 10^ 

T 



lOGH^^"^^ 



,5 \3/2 



(14) 
(15) 



For example, with lower resolving powers typical of spectrographs 
like VLT/UVES and Keck/HIRES (R = 70000), and assuming that 
(S /A/^)max = 300 can be easily obtained in a single comb exposure, 
a precision of cTy^^ = 2.4cms~^ should be achievable. The pre- 
factor of 0.45 cm s~^ in equations ([T4] ) & JTSl l is larger for narrower 
wavelength ranges. For example, it is 0.63 cm s~^ for 4500-6500 A. 

3.4 Additional sources of error 

The typical photon-limited precision of ~1 cm s~^ calculated above 
is 5 orders of magnitude smaller than the typical pixel size for most 
modern echelle spectrographs. To provide some context, for a 15- 
[im pixel this corresponds approximately to the size of the silicon 
atoms making up the CCD substrate! Therefore, if the photon noise 
limit is to be reached, mitigation of possible systematic eff'ects 
in the telescope, spectrograph and detector systems is imperative. 
Many such eff'ects have already been considered for observations 
requiring short-term stability (~ 1 night to ~ 1 week) which achieve 
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100-cms~^ precision or better (e.g. lLovis et al.ll2006h . However, in 
pushing to the l-cms~^ precision level, additional sources of error 
will need to be considered in detail. 

We will not conduct an exhaustive study of possible effects 
here but the following initial list is illustrative: (i) Inhomogeneities 
and variations in the intensity profile of the object and comb fi- 
bres presented at the spectrograph slit; (ii) Variations in the spectro- 
graph's instrumental profile (IP) across the CCD; (iii) Inaccuracies 
in the echelle grating's line markings; (iv) Variations in pixel size 
and separation across an individual CCD array; (v) Intra-pixel sen- 
sitivity variations; and (vi) CCD temperature variations, especially 
during read-out. While pre- slit effects such as (i) clearly require 
their own solutions (in this case, efficient beam homogenization 
and stabilization), most post- slit effects such as (ii-vi) might be 
mitigated by using the frequency comb itself. Since the comb pro- 
vides a series of closely spaced modes which are known, a priori, to 
have equidistant spacing, comb spectra could be used to accurately 
characterize the spectrograph and detector's response to different 
systematic effects. In practice, this would be facilitated by injecting 
the comb spectrum into the object fibre and by altering and/or 
Vce slightly to shift the comb spectrum by less than a resolution ele- 
ment. Thus, while the frequency comb would provide, in principle, 
a significant increase in precision, it would also provide an impor- 
tant means to characterize, measure and correct-for many system- 
atic errors so that the available photon-limited precision could be 
realised in practice. 

Aside from the spectrograph's IP, an additional factor con- 
tributing to asymmetries and variations in the measured point- 
spread function (PSF) could be intrinsic asymmetries in the comb 
modes themselves. Although the central frequencies of the modes 
are defined precisely by the self-referencing system, asymmetries 
in the modes' profiles could arise via correlated phase and am- 
plitude noise, generally referred to as 'repetition-rate noise'. The 
time-dependence (and therefore the time- averaging properties) of 
the asymmetry in individual modes will depend greatly on the type 
of laser comb system eventually employed. Nevertheless, since the 
spectral width of individual modes is typically ~1 MHz (0.6 ms~^ 
at 6000 A) and the asymmetry is expected to be small in comb 
systems stabilized to precise clocks, the difference between an in- 
dividual mode's frequency and the centroid of its intrinsic profile 
should amount to no more than ~l-10cms~^ If the asymmetries 
vary very slowly with time, this effect would be corrected for when 
modelling the overall PSF. However, if mode profiles were found 
to vary on very short time-scales, additional techniques for reduc- 
ing the mode line-width closer to the desired calibration uncer- 
tainty (i.e. -1 cms~^ or 17 kHz at 6000 A) might be employed (e.g. 



iDiddams et aDl2003h . 



4 DESIGN CHALLENGES AND POSSIBLE SOLUTIONS 

From the point of view of designing a real frequency comb cali- 
bration system, the most important result from the simulations in 
Section 13.31 is that laser repetition rates of Vr = 5-30 GHz give 
the highest velocity precision for practical spectrograph resolving 
powers ofR^ 50000-300000 [equation Although high repe- 
tition rate (> 10 GHz) lasers with picosecond pulses and frequency 
combs with broad spectral coverage, sub-lOO-fs pulses and low rep- 
etition rate (< 1 GHz) have been demonstrated, the desired com- 
bination of high repetition rates, short pulses and large spectral 
widths (e.g. ~3800-8200 A) is challenging and is yet to be experi- 
mentally demonstrated. A uniform intensity distribution across the 
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Figure 5. A review of femtosecond laser performance data shows two dis- 
tinct clusters, one in the 10-200-fs regime and the other in the 1-20-ps 
regime of pulse durations. Sub-lOO-fs pulse lasers seem to operate to up 
to 4 GHz. The region of interest (5-30 GHz repetition rate, 3-200 fs pulse 
duration) is still new territory for femtosecond pulsed lasers. 



whole wavelength range is also highly desirable since this ensures 
that similar velocity precision can be obtained from any spectral 
sub-division. Below we discuss why these design requirements may 
be problematic and offer some possible solutions. 



4.1 Repetition rate 

Reaching the optimum repetition rate suggested by Fig. |4] is not, 
by itself, challenging; several systems with even higher Vr already 
exist. For example jLecomte et al.l(l2005l) published a mode-locked 
laser with Vr = 40 GHz. Unfortunately, the peak power of such 
lasers is so low that they cannot use the fast but weak pulse shap- 
ing and shortening mechanisms that are usually provided by the 
Kerr effect. Instead these lasers rely on slow saturable absorbers for 
mode locking, forcing pulse durations into the picosecond regime 
with the associated narrow spectral widths. Indeed, the spectral 
width would need to be >30 times larger to maintain the comb's 
coherence in the spectral broadening process discussed below in 
Section |4l2l The challenge of generating sub- 100 fs pulses at repe- 
tition rates > 10 GHz is illustrated by a review of the performance 
of ultra-fast lasers in the literature, as summarised in Fig. (5] To our 
knowledge, the highest repetition rat e achieved so far in the sub- 
100-fs regime is v^- ^ 4 GHz (Leburn et al. 2004). 

One possible way to combine high repetition rate, short pulses 
and large spectral width is to spectrally filter the laser beam outside 
the cavity. The strategy would be to use a laser with modest repeti- 
tion rate, say Vr ~ 2 GHz, where reasonably short pulses can still be 
generated and where there is still enough energy per pulse to allow 
efficient spectral broadening (as discussed in Section |4l2l ). The ef- 
fective mode spacing could then be increased by filtering out most 
modes with a Fabry-Perot etalon with a free spectral range chosen 
to be an integer multiple of Vr and equal to the final desired mode 
spacing (e.g. lUdem et al .1 11999). The etalon would be required to 
have highly stable transmission fringes and a high enough finesse 
to ensure strong suppression of the extraneous modes. 

Figure [6] shows an etalon transmission function defined by the 
free spectral range (FSR), mode suppression ratio in decibels (pae) 
and by a repetition rate enhancement factor q = FSR/Vr. Since the 
cavity transmission profile should be adequately described by an 
Airy function, the required finesse of the etalon cavity is uniquely 
determined by pdB and q: 
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Figure 6. Simple representation of the (Fabry-Perot) transmission profile of 
an idealized mode-filter cavity (black curves). The resulting mode spacing 
is given by the free spectral range (FSR) between transmission peaks. The 
degree to which the intervening modes are suppressed is specified by the 
side-mode suppression ratio, pdB. measured at the position of the strongest 
ofF-resonant mode. Achieving strong side-mode suppression requires a high 
finesse according to equation {16}. 



2 sm(7T/q) 2 



(16) 



where the approximation is accurate for suppression ratios > 10 dB 
and ^ » 1. For example, the approximation is correct within 1.6 per 
cent for pdB = 50 dB and ^ = 10. For a 10-GHz FSR, a side-mode 
suppression of pdB = 50 dB can be achieved by eliminating 3 out of 
every 4 modes (i.e. q = 4) from a Vr ~ 2.5 GHz comb, requiring a fi- 
nesse of F ^ 700. Such a repetition rate and finesse are achievable 
with current technology and so this mode-filter technique seems 
a promising approach to achieve higher effective repetition rates. 
However, one challenge will be for a large spectral range to be fil- 
tered in this way because the dispersion of the cavity components 
will create differing path-lengths for light of widely different wave- 
lengths. The etalon mode- spacing remains constant only if the dis- 
persion is compensated properly. Simultaneously achieving reliable 
dispersion compensation and high finesse/reflectivity may impose a 
challenge for the etalon mirror design ( Gohle et al. 2005) . A prac- 
tical solution might employ several frequency combs and etalons 
which each cover a narrower spectral range. 



4.2 Spectral width 

In most octave- spanning frequency combs, the process of self- 
phase modulation is used to extend the comb's spectral width. For 
example, titanium-sapphire femtosecond laser pulses launched into 
photonic crystal fibres or micro- structured fibres have shown ex- 
tensive spectral broadening and this has been used for generat- 
ing frequency combs with Vr up to 1 GHz (e.g. Jones et al. 2000; 
iHolzwarth et alJbOOOl) . In this process the spectral coverage is in- 
creased fairly symmetrically around the central wavelength of the 
pulses and a total spectral coverage from 4500 A to 1.4iJ,m has 
been achieved. Similarly, lasers based on Er-doped fibres emitting 
around 1.5 ^tm have prod uced spectra covering 5000 A to 2.2 iJ,m 
at Vr up to 250 MHz fe.g. lNicholson et alllioosl : iMackowiak et all 
l2005h . 

However, producing broad spectra at higher repetition rates 
is difficult because the self-phase modulation process requires a 



certain pulse peak power to provide an octave of broadening. As- 
suming a pulse train at Vr with an average power, (P>, the peak 
power is roughly given by the energy in one pulse, (P) /Vr, divided 
by the pulse duration. Typically, a pulse energy of 2nJ for 100-fs 
pulses is needed, translating to average powers of 2 & 20 W for 1 
& 10 GHz repetition rates, respectively. In practice, (P) in mode- 
locked lasers have so far been limited to ~ 1 W. Thus, the high- 
est repetition rates achieved so far for octave- spanning combs are 
around Vr ~ 1 GHz. The fact that sufficient self-phase modulation 
requires a certain pulse peak power means that, in practice, broad 
spectra can be produced only with relatively low Vr while achieving 
higher Vr reduces the peak pulse power below that required for suffi- 
cient broadening . For ex ample, the repetition rate of 4 GHz demon- 
strated by Leburn et al.! (12004) is already close to optimal for very 
high spectrograph resolutions R ~ 300000 [equation ([T3])1. With a 
modest effort to filter the modes with an etalon, the effective repe- 
tition rate could probably be increased to > 10 GHz. However, this 
laser emits between 1.40 and 1.58 ^im so significant non-linear con- 
version processes would be required, such as frequency doubling 
(or tripling), or spectral shifting and broadening; the pulse peak 
power is too low at such high repetition rates for these processes to 
be effective. 

On the other hand, a new class of mode-locked lasers, opti- 
mized for larger spectral width rather than shorter pulses, is cur- 
rently under development. Such systems have reached spec tral 
widths extending from 5700 A to 1.3 ^im ( Miicke et al. ''2005") at 
repetition rates up to 1 GHz. However, even these broad spectra are 
not yet sufficient to cover the entire optical region, especially the 
blue region of the spectrum. The current limitation is actually due 
to the laser cavity mirrors which do not operate efficiently outside 
this (already relatively broad) range. 

Another option to achieve the desired effective spectral width 
could be the use of very fast and efficient optical phase modu- 
lators. These so-called 'optical frequency comb generators' (e.g . 
iKourod. Nakagawa & Ot"shulll993l : llmai Kourod & Otshulll998l) 
are capable of adding several hundred modulation side-bands to 
the spectrum of a single-mode laser with separations at the re- 
quired optimal mode spacing (i.e. ~ 10 GHz). However, so far the 
spectral width of these comb generators is insufficient for covering 
the whole optical range. Therefore, dozens (even hundreds) of gen- 
erators would have to be used simultaneously or some advanced 
scanning method would have to be devised which allowed the full 
calibration region to be swept with a single generator (or perhaps 
many) within the exposure time. 

In summary, designing a frequency comb covering the entire 
optical range while maintaining a high repetition rate would re- 
quire some improvements in comb technology. A practical solution 
might employ a combination of mode filtering stages and nonlinear 
conversion stages, possibly with additional amplifier stages to keep 
the power sufficiently high. 



4.3 Uniform intensity 

Unlike the sketch in Fig. [T] real frequency combs display strong 
intensity structure; intensity variations of 2 orders of magnitude 
are to be expected, even over relatively short wavelength ranges 
(~300A), especially when using photonic crystal fibres for spec- 
tral broadening (e.g. Holzwarth et al. 2001). The intensity structure 
is also sensitive to the polarization and peak power of the pulses 
coupled into the non-linear fibre. They therefore vary with time as 
the laser-fibre coupling and laser parameters change. Fortunately, 
since only a small amount of power is required to fully expose 
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the CCD (amounting to -40000 photons per CCD pixel over -10- 
3600-s exposure times), the spectrum can be flattened by an attenu- 
ation mechanism which can respond to both spectral and time vari- 
ations. One suc h mechanism , for example, might be a spatial light 
modulator (e.g. I Weine j|200d) whereby groups of modes (i.e. small 
chunks of the comb spectrum) can be separated, diff'erentially atten- 
uated and recombined again. An alternative approach is to improve 
the matching of phase velocities ('phase matching') between the 
various modes in the non-linear broadening process. This can be 
achieved by tailoring the longitudinal dispersion properties of the 
fibre and has been shown b oth the oretically and experimentally to 
produce broad, flat spectra (iHori et al. 2004) . 



5 CONCLUSION 

We have outlined the possible use of frequency combs as calibra- 
tion sources in astronomical echelle spectroscopy. Such combs gen- 
erate a series of equally spaced, very narrow modes which, by em- 
ploying self-referencing (or possibly other calibration) techniques, 
have absolute frequencies which are known a priori to relative pre- 
cisions better than 10~^^. These are significant advantages com- 
pared to the properties of other calibration sources such as ThAr 
emission-line lamps and h absorption cells. Simulations of comb 
spectra with uniform intensity over the optical range were sam- 
pled according to the properties of typical high resolution echelle 
spectrographs and CCD detectors. These revealed that the laser 
repetition rate must be Vr ~ 10-20 GHz to produce optimal line- 
spacings for R ~ 50000-200000 spectrographs. For such sys- 
tems, the photon-limited wavelength calibration precision is in the 
cms~^ regime when integrated over ~ 1000 A ranges and is gener- 
ally <10 cm s~^ when integrated over individual echelle orders. The 
regular grid of comb lines would also greatly facilitate the tracking 
and eff'ective removal of systematic distortions of the wavelength 
scale which might be produced through a variety of diff'erent eff'ects 
in the telescope, spectrograph and detector systems. The absolute 
wavelength calibration would therefore be stable over very long 
time- scales even if elements of the spectrograph or comb changed. 
Indeed, frequency comb calibration would place spectra observed 
on diff'erent telescopes at diff'erent times on the same absolute, high- 
precision wavelength scale. 

However, the results of the simulations reveal several chal- 
lenges in designing and implementing a working comb calibration 
system. Firstly, stable self -referenced combs with repetition rates of 
Vr ~ 10-20 GHz and which cover the entire optical range have not 
yet been realised. Another significant challenge will be to ensure 
that the comb envelope function - the peak intensity of the comb 
lines - is uniform over the wavelength range of interest. Never- 
theless, frequency comb technology is improving at a rapid pace 
and we expect that these challenges might be overcome by a rel- 
atively small research and development eff'ort. The potential gain 
is clear: frequency comb calibration may eff'ectively remove wave- 
length calibration uncertainties from all practical high-resolution 
spectroscopy and may allow the reliable combination of data taken 
on diff'erent telescopes over many decades. 

ACKNOWLEDGMENTS 

We are indebted to Savely Karshenboim for first suggesting the 
possibility of using frequency combs for wavelength calibration of 



astronomical spectra and to Robert Wynands for many early dis- 
cussions about the practicality of such a system. We thank Victor 
Flambaum and Jochen Liske for many informative and helpful dis- 
cussions. MTM thanks STFC for an Advanced Fellowship at the 
loA. 



REFERENCES 

Bahcall J. N., Sargent W. L. W., Schmidt M., 1967, ApJ, 149, Lll 
Baranne A. et al., 1996, A&AS, 119, 373 

Bauch A., 2003, Measurement Science and Technology, 14, 1159 
Bouchy P., Pepe P., Queloz D., 2001, A&A, 374, 733 
Brault J. W., 1987, Mikrochimica Acta, 3, 215 

Diddams S. A., Bartels A., Ramond T. M., Gates C. W., Bize S., Curtis 
P. A., Bergquist J. C, HoUberg L., 2003, IPPP J. Selected Topics Quan- 
tum Plectron., 9, 1072 

Portier T. M., Bartels A., Diddams S. A., 2006, Opt. Lett., 31, 1011 

Gohle C, Udem Th., Herrmann M., Rauschenberger J., Holzwarth R., 
Schuessler H. A., Krausz P., Hansch T. W., 2005, Nat, 436, 234 

Grazian A. et al. 2007, A&A, submitted 

Holzwarth R., Reichert J., Udem Th., Hansch T. W., 2001, Laser Phys., 11, 
1100 

Holzwarth R., Udem Th., Hansch T. W., Knight J. C, Wadsworth W. J., 

Russell R S. J., 2000, Phys. Rev. Lett., 85, 2264 
Hori T, Takayanagi J., Nishizawa N., Goto T, 2004, Opt. Express, 12, 317 
Imai K., Kourogi M., Otshu M., 1998, IPPP J. Quantum Plectron., 34, 54 
Ivanchik A., Petitjean P., Varshalovich D., Aracil B., Srianand R., Chand 

H., Ledoux C, Boisse R, 2005, A&A, 440, 45 
Jones D. J., Diddams S. A., Ranka J. K., Stentz A., Windeler R. S., Hall 

J. L., Cundiff S. T, 2000, Sci, 288, 635 
Knight J. C, Birks T. A., Russell R S. J., Atkin D. M., 1996, Opt. Lett., 21, 

1547 

Kourogi M., Nakagawa K., Otshu M., 1993, IPPP J. Quantum Plectron., 
29, 2693 

Leburn C. G., Lagatsky A. A., Brown C. T. A., Sibbett W., 2004, Plec- 
tron. Lett., 40, 805 
Lecomte S. et al, 2005, IPPP J. Quantum Plectron., 41, 45 
Loeb A., 1998, ApJ, 499, LI 11 
Lovis C. et al., 2006, Nat, 441, 305 

Mackowiak V., Kubina P., Adel P., Hansch T, Holzwarth R., 2005, in Quan- 
tum Plectron. Laser Sci. Conf. Vol. 2, QPLS'05. Opt. Soc. America, 
Washington, DC, U.S.A, p. 1188 

Marcy G. W., Butler R. R, 1992, PASP, 104, 270 

Marcy G. W., Butler R. R, 1996, ApJ, 464, L147 

Matos L., Kleppner D., Kuzucu O., Schibli T. R., Kim J., Ippen P. P., Kaert- 

ner P. X., 2004, Opt. Lett., 29, 1683 
Mayor M. et al., 2003, The Messenger, 1 14, 20 
Mayor M., Queloz D., 1995, Nat, 378, 355 

Mucke O. D., PU R., Winter A., Kim J.-W., Birge J. R., Matos L., Kartner 

P. X., 2005, Optics Pxpress, 13, 5163 
Murphy M. T, Plambaum V. V., Webb J. K., Dzuba V. V., Prochaska J. X., 

Wolfe A. M., 2004, Lecture Notes Phys., 648, 131 
Murphy M. T, Webb J. K., Plambaum V. V., 2003, MNRAS, 345, 609 
Murphy M. T, Webb J. K., Flambaum V. V., Dzuba V. A., Churchill C. W., 

Prochaska J. X., Barrow J. D., Wolfe A. M., 2001, MNRAS, 327, 1208 
Nicholson J. W. et al., 2003, Opt. Lett., 28, 643 

Pasquini L. et al., 2006, in Whitelock P., Dennefeld M., Leibundgut B., eds, 
lAU Symp. Ser. Vol. 232, Scientific Requirements for Pxtremely Large 
Telescopes. Cambridge Univ. Press, Cambridge, UK, p. 193 

Reichert J., Holzwarth R., Udem Th., Hansch T W., 1999, Opt. Communi- 
cations, 172, 59 

Reinhold P., Buning R., HoUenstein U., Ivanchik A., Petitjean P., Ubachs 

W., 2006, Phys. Rev. Lett., 96, 151101 
Rupprecht G. et al. 2004, in Moorwood A. P. M., lye M., eds, Proc. SPIP 

Vol. 5492, Ground-based Instrumentation for Astronomy, p. 148 
Sandage A., 1962, ApJ, 136, 319 



10 M. T. Murphy et al 



Udem Th., Holzwarth R., Hansch T. W., 2002, Nat, 416, 233 
Udem Th., Reichert J., Holzwarth R., Hansch T. W., 1999, Phys. Rev. Lett., 
82, 3568 

Varshalovich D. A., Levshakov S. A., 1993, J. Exp. Theor. Phys. Lett., 58, 
237 

Webb J. K., Flambaum V. V., Churchill C. W., Drinkwater M. J., Barrow 

J. D., 1999, Phys. Rev. Lett., 82, 884 
Weiner A. M., 2000, Rev. Sci. Instrum., 71, 1929 
Ye J., Ma L. S., Hall J. L., 2001, Phys. Rev. Lett., 87, 270801 



This paper has been typeset from a T^/ET^ file prepared by the author. 



